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Abstract. We report on an analysis of the drifting subpulses of PSR B0320+39 that indicates a sudden step of 
~ 180 degrees in subpulse phase near the centre of the pulse profile. The phase step, in combination with the 
attenuation of the periodic subpulse modulation at pulse longitudes near the step, suggests that the patterns 
arise from the addition of two superposed components of nearly opposite drift phase and differing longitudinal 
dependence. We argue that since there cannot be physical overlap of spark patterns on the polar cap, the drift 
components must be associated with a kind of "multiple imaging" of a single polar cap "carousel" spark pattern. 
One possibility is that the two components correspond to refracted rays originating from opposite sides of the 
polar cap. A second option associates the components with emission from two altitudes in the magnetosphere. 
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1. Introduction 

For a number of pulsars it is known that each pulse is 
composed of a number of subpulses, and that for each 
successive pulse the subpu lses appear to "drift" by a given 
amount across the profile l)Drake fc Crafdll968|) . The rate 
at which subpulses drift across the profile is not con- 
stant but rather depends on their pulse longitude. This 
can be understood in the context of the widely adopted 
"carousel" model of drifting subpulses ( Rudcrman 1972), 
which postulates the existence of a ring of equally spaced 
"sparks" above the polar cap that give rise to "tubes" of 
plasma streaming upward from the surface, ducted along 
magnetic field lines. At a certain altitude these particles 
emit microwave radiation that is beamed along tangents 
to the local field lines, giving rise to a beam pattern that 
consists of a ring of "subbeams" (bright spots), reflecting 
the polar cap spark configuration. As the pulsar rotates 
the distant observer samples emission along a certain path 
in the beam pattern. Under the right viewing geometry, 
the sight-line makes a tangential pass along the ring of 
subbeams, giving rise to the reception of a sequence of 
one to a few subpulses every time the star rotates. Due to 
the non-linear mapping between (spin) longitude and mag- 
netic azimuth, the subpulses should appear more closely 
spaced around the point at which the sight-line makes its 
closest approach to the magnetic pole. As the carousel 
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slowly rotates about the magnetic pole, the longitude at 
which the emission from a given spark is seen exhibits 
a monotonic drift with time, which by extension of the 
above argument is greater in magnitude further from the 
magnetic pole. 

By stacking the intensity time series modulo the pulse 
period to produce a two-dimensional array in pulse num- 
ber and pulse longitude, and taking Fourier transforms 
along constant-longitude columns (for ming a l ongitude - 
resolved fluctuation spectrum; LRFSl. iBacked l)l970albr ) 
was able to measure not only the characteristic peri- 
odicity in pulse number (P3, the time taken for one 
subpulse to drift to the former position of its neigh- 
bour), but also its amplitude and phase as a function 
of pulse longitude. Many subsequent studies employed 
this method but ignored the phase information, however 
those that examined the longitude-phase relation fBackei| 
1970a: Wright 1 981: Dayics ct al. 1984- .Biggs ct al. 1987; 
lAshworthingSSj) found behaviour qualitatively consistent 
with that expected under the carousel model. 

The longitude-phase relationship expected under the 
carousel model depends on the number of sparks present, 
the misalignment angle between the spin and magnetic 
axes, and their orientation with respect to the line of sight. 
Hence studies of the relationship have the potential to pro- 
vide valuable information about the pulsar, particularly in 
comparison to the observed longitudinal dependence of the 
polarisation position angle which, under the magnetic pole 
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model of lRadhakrishnan fc Cook3 l)l969l) . depends on the 
same geometric parameters. 

Motivated by these goals, we embarked on a program 
of analysis of several pulsars with drifting subpulses us- 
ing the Westerbork Synthesis Radio Telescope (WSRT). 
One of the target sources was PSR 603 20+39. This pulsar 
was discovered inde pendently by .Damashek et alJ l)l978|) 
and llzvekova et aD |i982) in surveys of the northern sky 
at 400 and 102.5 MHz respectively. Further observations 
revealed the em ission of very regular drifting subpulses 
lllzvekova et al ] 11982) which occur in two distinct pulse 
longitude intervals, wi th a non-drifting coniponent present 
in th e average profile l|Izvekoza et al, 1 11 993t 1 iKouwenhc 
l2000|) . In this paper we report a most unexpected result 
from our study of the longitude-phase relationship of the 
drifting subpulses of PSR B0320+39. 



2. Observations and Analysis 

We observed PSR B0320-f 39 on 2000 November 3 with the 
Westerbork Synthesis Radio Telescope (WSRT). Signals 
from fourteen telescopes were added (after appropriate de- 
lays) in each sense of linear polarisation in a 10 MHz band 
centred at 328 MHz. These were processed using the Dutch 
Pulsar Machine (PuMa), operating as a digital filterbank 
with 64 channels and 819.2 /is sampling in total inten- 
sity (S tokes I). For details of the system see lVoute et all 
ll2n02t) . 

In offline analysis we de-dispersed the data and 
binned the resultant time series into an array in pulse 
longitude and pul s e num ber using the cphemeris of 
lArzoumanian et alJ l)l994l) . A total of 3400 pulses in 3701 
longitude bins was recorded, with a region of 512 bins in- 
cluding the on-pulse used in further analysis. The samples 
were normalised to produce an average profile with a zero 
baseline and a peak value of unity. The first 100 pulses are 
shown in a greyscale form in Fig. 

The techniques employed in the further analys is are 
based on those presented in detail by^dwards & Sta ppergj 
ll2002l) . We computed the two-dimensional Discrete 
Fourier Transform (DFT) of the first 2048 pulses of the 
pulse longitude - pulse number array to yield what we call 
the two-dimensional fluctuation spectrum (2DFS). This 
is mathematically equivalen t to the harmonic resolved 
fluctuation spectrum (HRF S; iDeshpande k. Rankinlll999t 
lEdwards fc Stappersll200^ ^ and the form of the power 
spectrum found here fo r PSR B0320-I-39 matches that 
found by iKouwenhovenI ((200nl) for the HRFS. For the 
purposes of clear display we also computed a spectrum 
using data from which we subtracted the average pro- 
file from each pulse, and padded with zeros in the lon- 



Without the use of "twiddle factors" dCoolev fc TukevI 
[l96^, the 2DFS will differ from the HRFS. However, it is eas- 
ily shown that this merely amounts to the provision of slightly 
different sampling lattices over the same continuous function, 
due to the assignment of integer pulse numbers (2DFS) versus 
integer harmonic numbers (HRFS). 
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Fig. 1. First 100 pulses of 328 MHz timeseries for PSR 
B0320+39, binned in pulse longitude and pulse number. 



gitude axis, in order to improve the resolution. The reso- 
lution in the time-associated axis was unnecessarily high 
(2048 elements), so to improve the signal-to- noise ratio 
we smoothed the power spectrum by convolution with a 
7x1 -bin {vt x vi) boxcar before plotting the relevant 
portion of the spectrum^ in Fig. (21 The spectrum is dom- 
inated by a pair of components associated with the drift 
modulation; two are present because the spectrum de- 
rives from the Fourier transform of real-valued data (see 
lEdwards fc Stappers,,2002i) . Also present with high signif- 
icance are a pair of responses corresponding to the sec- 
ond harmonic of the subpulse response. The features near 
Vt = 0.02 and 0.04 cycles/period are significant, however 
they also appear in an analysis of off-pulse data, indicating 
that they are caused by periodic interference. 

After performing the 2D DFT, the complex spectrum 
was shifted in both axes by an amount which placed the 
primary response to the drifting subpulses at zero fre- 
quency. We then masked other components present in the 
spectrum (namely the complex conjugate "mirror" of the 
primary response, the pair of second harmonics and the 
DC component) by multiplication of the coefficients with 
a transfer function: 




\vt - VtcMw < 0.5 
vtc\/w - 1 0.5 < \iyt - vtcMw < 1 , (1) 

\vt - Vtc\/W > 1 



where vt is the time-associated axis, vtc is the centre fre- 
quency of the component to mask, and w is the frequency 
extent of the region of zero transmission. The results that 



^ The sampling interval used allows longitude-associated fre- 
quencies up to 1850 cycles/period, however all significant 
power is present in the portion of the spectrum shown. 
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Fig. 2. Main panel: Two-Dimensional Fluctuation Power 
Spectrum of PSR B0320+39. The greyscale saturates at 
0.05 times the peak power, while the contours show the 
same spectrum at levels 0.05 and 0.5 times the peak. Left 
panel: "collapsed" spectrum formed by averaging rows of 
the 2D spectrum. Bottom panel: power as a function of vi 
for i/t = -0.118. 

follow derive from filters with w = 0.05 cycles/period, al- 
though as expected they were found to be insensitive to 
w unless it becomes too small to mask adequately or so 
large that it masks some of the desired signal. We then 
computed the inverse DFT of the result. 

We represent the subpulse signal as the real part of 
the product of a pure two-dimensional sinusoid and a two- 
dimensional complex modulation envelope. By the convo- 
lution theorem, this envelope is the result of the inverse 
transform performed above. It describes the deviations 
from a pure periodicity of constant amplitude that arise 
due to amplitude modulation at the pulse period, scin- 
tillation or nulling, non-uniform spacing of subpulses in 
pulse longitude (e.g. due to sight-line curvature) and slow 
variatio ns in the drift rate. Using the iterative scheme 
of Edwards fc Stapperj 12002) we decomposed the two- 
dimensional envelope into the product of two modulation 
envelopes that vary in pulse longitude and pulse number 
respectively. Such a decomposition is able to model the be- 
haviour of a rotating carousel, where the geometry defines 
the longitude-phase dependence and the carousel rotation 
determines the time-phase dependence, and its basic as- 



sumptions - that the longitudinal spacing of the subpulses 
is time invariant, and that the time spacing of the sub- 
pulses is longitude i ndependent - are in agreein ent with 
previous studies; see lEdwards fc Staopersl (jioQ^)- The va- 
lidity of this approach in the case of PSR B0320-I-39 was 
confirmed by computation of the 2DFS of the difference 
between the observed signal and that predicted using the 
real part of product of the two one-dimensional envelopes. 
No significant power remained in the portion of the spec- 
trum previously occupied by the primary drifting subpulse 
response. 

The spectral shifting performed earlier has the effect of 
removing a constant phase slope in each of the envelopes. 
Since the slope varies over the envelope, the choice of value 
to use in the shifting is somewhat arbitrary and was re- 
fined after viewing the resultant envelopes. The period- 
icity with pulse number was well defined and a value of 
P3 — — Pi/0.118 was used. The response to the drifting 
component in the other axis is very broad, and in fact 
extends well into negative frequencies (Fig.[21). This com- 
plicates the choice of value for I/P2 to be used in shifting 
the spectrum in the longitude axis. We chose a value of 
P2 = Pi/QO to give the resultant envelope a relatively 
small phase slope over most of the pulse^. 

The inferred longitude-dependent envelope (Fig. |2Jl 
shows that nearly linear subpulse modulation occurs in 
two distinct longitude intervals, with a striking phase off- 
set of ^180° between them. In terms of the longitude-time 
diagram, the effect of the offset is to cause the drift band 
peaks in one component to extrapolate to troughs in the 
other and vice versa. To check this result we computed 
the (complex) LRFS and examined the amplitudes and 
phases in frequency bins near P3 — Pi/0.118. Since the 
modulation is not a pure periodicity, the power is spread 
over several bins in the spectrum. Although of a much 
lower signal-to-noise ratio, the shape of the phase enve- 
lope in each bin was consistent with that inferred with our 
preferred technique, which makes optimal use of all mod- 
ulation power. Similarly, division of the observation into 
segments in time (which results in wider frequency bins 
and potentially also in reduced frequency drift in each seg- 
ment) and examination of the phase relation in the peak 
bin of each segment also gave consistent results (regard- 
less of the segment size). In order to improve the signal- 
to-noise ratio of these other methods, we also used an iter- 
ative algorithm based on multiplication with the complex 
conjugate of a template envelope to determine the correct 
phases with which to add envelopes from different fre- 
quency bins or time segments (producing in effect coa rse- 
grain ed versions of the algorithm of Edwards & Sta pperj 
I2OO2I) . The resultant summed envelopes were again con- 
sistent with the result of our preferred algorithm (Fig. ISJ . 
As a final check, we folded the longitude-time data about 
the pulse number axis (modulo 1/0.118 pulses) to produce 

^ Since this value does not lie on the natural grid of the DFT, 
the shift was actually performed in the time-longitude domain 
by multiplication with a complex exponential. 
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Fig. 3. Inferred modulation envelope (solid line and filled 
circles) and average pulse profile (dotted line). Also shown 
is the phase envelope computed by coherent addition of a 
set of 69 rows (0.108 < vt < 0.128) of the 3400-pt LRFS 
over which significant power was seen (unfilled circles), 
and the phase envelope computed by coherent addition of 
the peak row from each of twenty-six 128-pt LRF spectra 
resulting from time-segmentation of the sequence (with 
this transform length the majority of modulation power 
falls in a single row). The nominal phase slope of I/P2 — 
60°/° has been subtracted from the phase envelopes before 
plotting. 



an average profile as a function of drift phase and pulse 
longitude (Fig. 0J. As with the LRFS methods, much of 
the modulation was lost due to the unstable drift motion, 
however some residual modulation remained and is suffi- 
cient to confirm the unusual phase behaviour. 

Knowing that there is a ~180° offset between the two 
halves of the subpulse phase envelope, we may explain 
two otherwise curious facts about PSR B0320-I-39. The 
first concerns the val ues of P2 previously r eported, which 
are around Pi/120 iflzvekoza et al.l flQQSl report 3.1° of 
longitude, i.e. P2 — Pi/116 at 406 MHz; iKouwenhovenI 
[2000, reports 24.5 ms, i.e. P2 ~ 123, at 328 MHz). These 




Pulse longitude (deg) 

Fig. 4. Data folded at P3, as a function of pulse longitude 
and drift phase. Maxima and minima in each longitude bin 
are indicated with black and white crosses respectively. 
The data are plotted four times to guide the eye. 



strongly differ from our nominal P2 = Pi /60 because they 
arc based on measurements of actual subpulse separations, 
which will nearly always be taken between subpulses on 
opposite sides of the phase shift. Since the subpulse in the 
trailing component leads (or lags) the position it would 
take in a continuous linear phase envelope by ~ 180° of 
subpulse phase, the derived spacing differs from that ex- 
pected from the typical phase slope by a factor of 1 /2 (and 
probably sometimes 3/2). The second curiosity explained 
by the phase offset is the local minimum in the subpulse 
HRFS response at around 60 cycles/period (Fig.[2Jl. The 
complex longitude-dependent modulation envelope (once 
the constant phase slope is subtracted) is quite close to an 
odd real-valued function, which causes its Fourier trans- 
form to have local minimum in its power at DC. This is 
the cause of the local minimum around 60 cycles/period, 
the typical longitudinal drift frequency. Without the phase 
offset the envelope would be close to a real-valued even 
function and its Fourier transform would not have a strong 
minimum at its centre. 

Furthermore, we note that the equivalence of the 2DFS 
and HRFS provides an explanation for the apparent "sec- 
ondary" component at 0.118 cycles/period, suggested by 
IKouwenhovenI l|200Cl() to be due to the ambiguity between 
which drift band joins to which between the two com- 
ponents. In fact, the subpulse response is very broad in 
the longitude- associated frequency axis, and extends sig- 
nificantly into negative frequencies. The symmetry of the 
two dimensional Fourier transform of real functions is such 
that the power at any frequency pair i^i , i^t is equal to 
that at —vi,—vt (where vi and vt are the frequencies in 
the longitude-associated and time-associated axes respec- 
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Fig. 5. Inferred subpulse phase envelope, including the 
60°/° phase slope previously removed (indicated by the 
straight line). The data are plotted twice with an offset 
of 180° to facilate evaluation of the central feature as the 
transition between two offset regions of linear drift. Two- 
sigma error bars are included. The curved line represents 
the best fit to the data under the standard carousel model, 
given the polarimetric constraint (see text). 

tively). If only the positive frequencies in vi are plotted, 
as is the case with the HRFS, no information is lost but 
individual components that span the i^; = axis will ap- 
pear as two spatially separated components at vt values 
of opposite sign. When the full spectrum is plotted as in 
Fig. 121 the proper interpretation becomes apparent. 

3. Discussion & Interpretation 

3.1. Polar Cap Models 

The phenomenon of drifting subpulses has been used re- 
peatedly in investigations of pulsar emission as a po- 
tentially strong clue about the geometry of the system. 
The fact that those pulsars that show nearly linear drift- 
ing subpulses tend to have a double, "unresolved dou- 
ble" or "single" average pulse profile morphologies, with 
many showing a transition from resolved to unresolved 
double as the observing frequency is increased, is usu- 
ally seen as strong supp ort for some form of "hollow- 
cone" polar cap model llRadha krishna n fc 



a transition from a resolved double to an unresolved 
double average profi le morphology between 102 MHz 
lllzvekova et all jl982l) ' | and 328 MHz (th is work; also 
iDamashek et a1.lll 97SHTzvekoza et h}\[^99^. 

The local steepening of the phase envelope of PSR 
B0320-f39 in the transition region is inconsistent with the 
carousel model in its most basic form. For a single ring 
of "sparks" , as long as the sparks are elongated in either 
magnetic azimuth or opening angle (or neither, but not 
both) , the subpulse modulation seen in a ring of constant 
opening angle has the same phase as that in any other ring, 
if measured from the same fiducial azimuth. Therefore, the 
observed longitude-dependent subpulse phase envelope is 
directly tied to to the magnetic azimuth sampling of the 
sight-line (via multiplication by the number of sparks, N) 
((Edwards fc StaDDersll200 j) . The only way to produce a 
sudden swing in magnetic azimuth (and hence subpulse 
phase) is to pass close to the magnetic pole. Such a con- 
figuration is inconsistent with the usual interpretation of 
unresolved double profile morphologies, and with the grad- 
ual linear polarisation position angle swin g observed in 
this pulsar i Suleimanova fc Pugachevl20n2l) . To prove this 
assertion we attempted to fit the observed subpulse phase 
swing with the geometric formula of Edw ards & Stapperty 
l|2002() , usin g: a minimisation algorithm. Fro m the figure 
provided bv ISuleimanova fc PugachevI l|2002() we conser- 
vatively estimate the range of polarization position angle 
gradients consistent with the data as dx/d(f) = 12 ± 3°/° 
(where x and (j) are the position angle and pulse longitude 
respectively). The parameters of the fit were required to 
be consistent with this grad i ent vi a the geometric model 
of iRadhakris hnan fc Cook3 l|l969l) . The result (Fig. 
clearly fails to describe the observed behaviour. 

One way to produce a jump in subpulse phase is 
to adopt a nested ring polar cap spark configuration 
('Oil fc Sendvk' '2nnm and to assign the two halves of the 
profile to rings of different opening angles and numbers of 
sparks. However, this would result in a measurably differ- 
ent P3 unless the circulation rates of the two carousels 
were different by a fortuitous ratio (i.e. the inverse of 
the ratio of the numbers of sparks in the rings). In ad- 
dition, the transition from resolved to unresolved double 
profile with increasing frequency could not be explained 
by st andard radius-to-frequency mapping (REM; ICordesI 
'l97^. Some other mechanism such as the finite beam- 
width of the elementary emission mechanism (27"-'^), the 
finite band-width of emission at a given radius (with at- 
tendant REM in the effective range of radii visible), or 
refractive broadening would be required to produce the 
profile evolution. 



i Ruderman fc Sutherlandll975HBackeill976HRankinll983 , 
Il98fit iLvne fc Manchester! Il98al) . PSR, B0320-(-39 shares 
these characteristics of regular drifting subpulses and 



^ Due to an imp r oved system response the profile of 
iKuz'min fc Loso vskii ( 1999) shows much stronger component 
separation than that of .lzvekova et aL (.1982.') : see fig-IHl 
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3.2. Double Imaging of Magnetospheric Origin 
3.2.1. Basic Considerations 

Given the failure of the standard model to account for the 
observed subpulse drift behaviour of PSR B0320+39, we 
develop in this section a pair of explanations based on a 
simple and plausible extension of the model. Namely, we 
argue that the most important features of the longitude- 
dependent subpulse modulation envelope - a local re- 
duction of amplitude accompanied by a rapid swing in 
the phase angle - are the result of destructive interfer- 
ence between two superposed drifting subpulse signals of 
nearly opposite phase. The observed behaviour and this 
explanation of it can be understood by analogy with the 
phenomenon of orthogonal polarisation mode transitions, 
where the pulsar signal is seen to consist of the incoher- 
ent sum of two signals with nearly orthogonal polarisa- 
tion states (e.g. [Manche ster et aLlllOTfil IStinebring et alJ 
In both cases the total intensity of the two signals 
simply add in the average profile, while the near anti- 
parallel nature of the vectors representing the complex 
subpulse modulation signal (in the Argand plane) or the 
polarised part of the signal (in the Poincare sphere) ex- 
plains the reduction in magnitude (subpulse amplitude, 
polarised intensity) and the swing in orientation (subpulse 
phase, polarisation position angle) during a transition of 
dominance from one mode to the other. Other than the 
obvious features of the shift in phase and the attenuation 
in amplitude, the continuity of both the average profile 
and the absolute value of the gradient of the subpulse 
amplitude envelope serve as supporting evidence for this 
hypothesis. 

Since the polar cap cannot be physically populated 
with overlapping spark systems, we argue that the super- 
posed drifting components of PSR B0320-I-39 are associ- 
ated with a form of "double imaging" of the polar cap 
excitation pattern. Whereas the polar cap may consist of 
a single, drifting ring of excitation points (e.g. sparks), the 
physics that map this to a radiated beam pattern must do 
so under two different transformations that are received in 
incoherent addition by the observer. This gives rise to a 
subpulse pattern that is the sum of two overlapping com- 
ponents with differing peak pulse longitudes, and phases 
that differ in their contribution at a given pulse longitude 
by a constant offset of ~ 180°. This also implies that the 
average pulse profile is the simple sum of two overlap- 
ping components approximately equal to the correspond- 
ing subpulse amplitude envelopes. 

In this section we discuss two potential mecha- 
nisms for such an effect. The first associates the images 
with different rav paths iii a refractive magnetosphere 
l|Petrova & LvubarskiilEKKil) . while the second attributes 
them to emission originating at two discrete altitudes in 
the magnetosphere (jRankln 1993) . Before going into the 
details of each scenario, given that they both produce 
a beam consisting of a pair of nested cones, we are are 
in a position to examine the consistency between these 
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Fig. 6. Average total intensity pulse profiles of PSR 
B0320-f 39 over a frequency range of 102-1408 MHz. The 
profiles are normalised to the same peak intensity and (in 
the absence of timing information) were aligned by eye 
to match the model proposed here. The 328 MHz pro- 
file is from this wor k . The 102 MHz profile is that of 
iKiiz' min & Losovskii ('1999'), and the remaining profiles 
are from Gould & Lvnc ( 1998). All profiles other than the 
328 MHz one were obtained via the EPN profile database. 

models and the observed pulse profile morphology. This 
is important because, in assigning the two components of 
PSR B0320+39 to cones of different opening angles, we 
are abandoning the standard "single cone" explanation 
favoured for most pulsars with regular drifting subpulses. 
We suggest instead that a cut through an outer cone (from 
the outside to the inside) produces a leading component, 
while a trailing component is produced by a pass along 
the edge of the inner cone. At most radio frequencies the 
difference in the opening angles of the cones is small and 
the components overlap, but radius-to-frequency mapping 
causes the outer cone to expand at lower frequencies, giv- 
ing the components sufficient separation to be resolved at 
100 MHz. Clearly one should also expect a third compo- 
nent to be present on the trailing edge of the profile as 
the sight-line passes across the outer cone for the second 
time. This is actually observed at frequencies at and above 
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600 MHz l|Gould fc Lvnelll998|) . giving further support to 
our assertion that more than one cone of emission is visi- 
ble from this pulsar. As shown in Fig|^ the pulse profiles 
observed for PSR B0320+39 match this picture well. We 
now turn to the specifics of the two mechanisms. 



nner cone 
ray 



3.2.2. Refraction 

Refractive effects have long been conside red poteiitially 
important in pulsar ni agnetospheres (e.g. iMelrosd Il979t 
iBarnard fc Arondl 19861). particularly in the p roduction of 
circular polarisation ijAllen fc Melrosd 119821 ) and super- 
posed orthogonally polarised radiation (McKinnon 1992). 
Only recently, however, has the potential for refraction 
to explain pulse profile morphologies seen detailed study 
('Lvubarskii fc Petrovalll998t IPetrova fc Lvubarskiill200nt 
[Petrova.2000(l . Based on the simple argument that the 
magnetospheric plasma density has minima both at the 
magnetic pole and at the edges o f the cone defined by 
the la st open magnetic field lines, IPetrova fc Lvubarskiil 
l)2000|) showed that rays emitted poleward of the peak 
of the plasma distribution could be refracted across the 
pole to exit with an angle to the pole of the opposite sign 
to that with which they were emitted. Rays emitted out- 
ward of the plasma peak are simply refracted to somewhat 
greater opening angles. This is illustrated schematically in 
Fig. 13 Clearly, at any pulse longitude the observer may 
potentially receive two rays, one from either side of the 
magnetic pole. 

In an axisymmetric plasma distribution, refra ction oc- 
curs p urely in the plane of the magnetic field lines. IPetroval 
l)200 0l) shows that the magnetosphere tends to concentrate 
rays from within a given set of co-planar field lines into two 
components of different peak opening angle, correspond- 
ing to the different ray paths described above. These are 
associated with inner and outer components of the conal 
beam, which in terms of explaining the morphology of pul- 
sar average profiles serves the same purpose as the dou- 
ble cone of lRankinI l)l993|) (see below). This grouping of 
rays occurs within a given plane of field lines; all possi- 
ble origins for visible rays must be co-planar, along with 
the line of sight and the magnetic axis. Therefore, to the 
extent that the plasma distribution is axisymmetric, the 
geometric conveniences of the b asic magnetic pole model 
( Radhakrishnan & Cook§ Il969() are preserved, including 
the polarisation behaviour, and the mapping of azimuthal 
intensity modulations (e.g. from a plasma flow initiated 
by a polar cap spark carousel) to observable drifting sub- 
pulses. 

Moreover, if the polar cap pattern consists of a ring of 
an odd number of points of excitation, the subpulse pat- 
tern seen in the component associated with pole-crossing 
rays should be close to 180° out of phase with that in 
the other component. The slight deviation from purely 
out of phase signals required to explain the finite longi- 
tude interval over which the phase transition occurs might 
be explained as a retardation effect or a refractive de- 




Fig. 7. Schematic illustration of the simultaneous observ- 
ability of rays from opposite sides of the polar cap, via 
magnetospheric refraction. The plane of the page contains 
the magnetic axis, the observer, and all possible origins for 
observable rays emitted and refracted under axisymmetric 
conditions. The "active" field lines that mark the peak of 
the plasma distribution at a given altitude are labelled. 
The field lines defining the approximate edge of the emit- 
ting flux tubes are shown with dotted lines. Two rays, 
corresponding to those seen by an observer whose sight- 
line presently intersects the overlap region of two conal 
beams, are shown. The left hand ray originates outside 
the peak of the plasma distribution, and is refracted fur- 
ther in this direction. The right hand ray originates inward 
of the peak plasma density and is refracted across the pole 
to exit at the same angle as the other ray. A carousel-like 
intensity pattern with 5 subbeams is shown to illustrate 
the origin of the subpulse phase difference between the 
two components associated with these ray paths. 
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Fig. 8. Schematic illustration of the simultaneous observ- 
ability of rays from two different emission heights (see also 
Fig. 0). Due to the curvature of the magnetic field lines 
along which the radiating particles stream, the ray from 
the upper region must come from a more inward field line 
that that from the lower region, if the two are to make the 
same angle to the magnetic axis. 



lay. Alternatively, inclusion of magnetospheric rotation ef- 
fects or refraction un der a moderately non-axissym metric 
plasma distribution IPetro va fc Lvubarskiil l2000l model 
both but only in t he context of polarisati on wave mode 
coupling; see also iBlaskiewicz et al.l Il99ll) could poten- 
tially provide the kind of symmetry-breaking needed to 
produce subpulse phase offsets other than or 180°. The 
form of the beam pattern and the manner in which a rotat- 
ing carousel polar cap system maps to drifting subbeams 
under such conditions is yet to be investigated, to our 
knowledge. 

We note that this model can also explain a recent re- 
sult concerning the drifting subpulses of PSR B0809-I-74. 



iRamachandran' et all l)2002j) found that over the course of 
each subpulse, the polarisation state jumps in a regular 
manner from one polarisation state to another, almost or- 
thogonal state. This is similar to the behaviour seen in 
PSR B0320+39, in that the observed total intensity drift- 
ing subpulse pattern can be understood as the sum of two 
components with a phase offset between them. However, in 
PSR B0320-I-39 the drift components cannot be associated 
with orthogonal polarisation modes since the linearly po- 
larised average profi le at 408 MHz peaks i n intensity at the 
centre of the profile llCxOuld fc LvnJl998^ ■ where maximal 
depolarisation ought to take place, were the components 
or thogonally polarised. The single- pulse polarimetric work 
of lSuleimanova fc Piigach'evI l|2002l) confirms that the situ- 
ation is similar at 103 MHz: significant power is p resent in 
two modes only in the profile wings. The model of lPetroval 
I2OOI) for the origin of orthogonal polarisation modes of- 
fers a simple explanation for this disparity. Under that 
model the secondary polarisation mode is produced via 
the conversion of power from the ordinary mode in regions 
of propagation that are quasi- longitudinal to the magnetic 
field. If radiation in one of the conal components experi- 
ences significant conversion while the other does not, as 
could well occur since the ray trajectories involved are 
very different between the two, the situation is as seen in 
PSR B0809-I-74. If the configuration is not conducive to 
conversion then both components have similar polarisa- 
tion, as in PSR B0320-I-39. Since the conversion process 
also alters the orientation (i.e. position angle, ellipticity) of 
the modes, the longitude-dependent non-orthogonality of 
the po larisation modes reported bv IRamachandran et alJ 
l|2002|) might also receive an explanation via this effect. 

3.2.3. Dual Emission Altitudes 

A second potential explanation for how we receive the sum 
of two related subpulse signals is that the two signals arise 
at different altitudes in the magneto sphere. This is the 
essence of the "double conal" model of lRankinI l)l993j) . so- 
called because the production of radiation cones at two 
distinct altitudes in the magnetosphere gives rise to a 
beam pattern that consists of two nested cones. The radial 
separation of the emitting regions and the differential tan- 
gent angle between the active field lines and the magnetic 
axis that arises as a result causes the polar cap plasma 
flux cone to produce a wider radiation beam at the up- 
per emission region, compared to the lower one. If the 
height difference is large, up to four well-resolved conal 
profile components are produced, while reduced vertical 
separation causes the components to overlap and produce 
a "boxy" profile shape. Fig. |H1 illustrates this model as 
it might apply for PSR B0320-t-39 at 102 MHz, with the 
conal beams overlapping to a lesser degree than at higher 
frequencies. 

Clearly this model is able to explain how we see su- 
perposed drifting patterns, but we have not yet offered 
an explanation for why there should be a subpulse phase 
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offset between them. They arise from the same pattern 
of sparks (or plasma flux production) on the polar cap, 
yet we observe a phase difference of ^ 180° near the 
profile centre. For the subpulse pattern to be delayed 
by P 9./2 ~ 30 ms due to simple retardation l|Cordej 
Il978l) requires that the emission regions are separated 
by ~ 9000 km, 1-2 orders of magnitude greater than 
the yalues typically inferred frq r n profile morphologies 
fe.g. ICangadhara fc CuntalEiTl iMitra, T^ankinl Eioi 
ICunta fc Gangadharall2003 ) . Including differenti al aber- 
ratio n reduces this by no more than a factor of 2 ijCorded 
Il978|) . Hence we consider this explanation unlikely. An al- 
ternatiye explanation worth pursuing in the future, but 
outside the scope of this work, is that a significant electric 
field exists throughout the magnetosphere, causing par- 
ticles to accelerate upward mostly along magne tic field 
lines, but with some azimuthal E x B drift (e.g. IWrightl 
1200311 . The drift would cause the beam pattern radiated 
at a giyen height to be rotated in magnetic azimuth by a 
certain amount, which may account for the subpulse phase 
offset seen here. 

In any case, under this model the phase difference must 
ultimately be due to the different altitudes at which the 
components originate. Since RFM must apply to the outer 
(i.e. upper) component to explain the 102 MHz pulse pro- 
file, the difference in subpulse phase between the compo- 
nents ought to be a strong function of frequency. Although 
this is the first study to examine the longitude depen- 
dence of subpulse phase in this pulsar, the fact that the 
longitude-dependent subpulse amp litude enyelope s mea,- 



sured at 102 MHz and 408 MHz bv llzyekoza et all ((1 
drop to close to zero in the profile centre, just as at 
328 MHz, suggests that at all three frequencies the sub- 
pulses add destructively in this interval, and therefore that 
the phase difference between them is close to 180°, regard- 
less of frequency. The origin of the phase offset therefore 
remains a mystery under this scenario. 

4. Conclusions 

The drifting subpulses of PSR B0320-I-39 occur in two 
distinct intervals of pulse longitude. Using a sensitive new 
technique we have shown that the extrapolation of sub- 
pulse drift from a given longitude interval is ~ 180° out of 
phase with the subpulses of the other. This was confirmed 
by analysis of the longitude-resolved fiuctuation spectrum 
and by folding the data into profiles corresponding to dif- 
ferent phases of the subpulse drift. 

This suggests the presence of two components of nearly 
opposite drift phase, and based on the forms of the aver- 
age pulse profile and the longitudinal dependence of the 
amplitude of subpulses, we argued that the modes are su- 
perposed at all longitudes. The adding of intensity contri- 
butions from the two modes results in attenuation of the 
modulations in the longitude interval where the modes 
are of similar intensity, and produces a steep transition of 
~ 180° in the longitudinal subpulse phase envelope at the 
point where the component intensities are equal. 



Two specific models are suggested to explain this be- 
haviour. Both produce a nested pair of conal beams, one 
via magnetospheric refraction, the other via emission at 
two distinct altitudes. The leading and trailing compo- 
nents are associated with outer and inner cones respec- 
tively. The second cut of the outer cone does not give 
rise to a component at 102 MHz or 328 MHz, but does 
produce a third, trailing component at frequencies above 
600 MHz. The component separation increases rapidly at 
frequencies below 328 MHz, in accordance with the ex- 
pectations of radius to frequency mapping. The refractive 
scenario can explain the observed phase offset as a simple 
consequence of having an odd number of subbeams in the 
carousel, whilst under the dual-altitude emission model 
the phase offset is more difficult to understand. Planned 
observations of the time and frequency dependence of the 
phase offset, and the behaviour of subpulses in the third 
component present at and above 600 MHz, should prove 
most helpful in the further evaluation of these models. 
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